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Abstract. Results are presented from a multi-wavelength study of the giant pillars within the Carina Nebula. Using 
near-IR data from 2MASS, mid-IR data from MSX, 843 MHz radio continuum maps from the MOST and molecular 
line and continuum observations from the SEST, we investigate the nature of the pillars and search for evidence 
of ongoing star formation within them. Photodissociation regions (PDRs) exist across the whole nebula and trace 
the giant pillars, as well as many ridges, filaments and condensations (A v > 7 mag). Morphological similarities 
between emission features at 21 /^m and 843 MHz adjacent to the PDRs, suggests that the molecular material 
has been carved by the intense stellar winds and UV radiation from the nearby massive stars. In addition, star 
forming cores are found at the tips of several of the pillars. Using a stellar density distribution, several candidate 
embedded clusters are also found. One is clearly seen in the 2MASS images and is located within a dense core 
(G287.84-0.82). A search for massive young stellar objects and compact Hn regions using mid-IR colour criteria, 
reveal twelve candidates across the complex. Grey-body fits to SEDs for four of these objects are suggestive of 
OB-stars. We find that massive star formation in the Carina Nebula is occurring across the whole complex and 
confirm it has been continuous over the past 3 Myrs. 

Key words. ISM: structure; ISM: lines and bands; ISM: molecules; Hn regions; dust, extinction; Stars: formation 



1. Introduction 

In many cases, young stars have been found at the tips 
• of giant pillars that point toward a more evolved mas- 
1 sive sta r cluster (e.g. the elephant trun ks of the Eagle 
. Nebula: IMcCaughrean fe Andersen! l2002t) . The formation 
' of such pillars can readily occur if a dense core within a 
giant molecular cloud (GMC) is exposed to the intense 
stellar winds and radiation fields from a nearby massive 
star cluster. The core would shield the column of molecu- 
lar material behind it, in a direction pointing away from 
the cluster. Subsequently, the more exposed parts of the 
GMC would be swept up around this column or be com- 
pletely irradiated away. 

It is not clear if such a drastic change in the structure 
of a GMC can affect its star formation capacity. There 
is growing evidence to suggest that the tips of pillars are 
prime sites fo r ongo i ng star formation l|jiang et aTir2 002: 
IStanke et ail l2002t IMcCaughrean fc Andersen! l2002|h 
However, there is much debate over whether this type of 
star formation has been triggered by external processes, 
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or whether it has spontaneously formed. It is uncertain if 
we can even distinguish between the two. 

Several giant pillars pointing toward the massive clus- 
ters within the Carina Nebula have recently been dis- 
covered at mid-inf rared (mid-IR) wa velengths, the largest 
extending ~25pc Ipmith et alJ 2000). Bright IR emission 
condensations are located at the tips of several of these 
pillars, and may correspond to sites where star formation 
has been triggered by the intera ctions with the sur round- 
ing young massive star clusters ijSmith et al.ll200f)(h 

Many large molecular clouds are associated with the 
Carina Nebula llZhang et al l 120011 iBrooks et all Il99 
IWhiteoak fc Otruncekl Il984t Ide Graauw et all Il98 
These clouds lie close to several massive star clusters, in 
particular Bochum 10 and 11, Collinder 228 and Trumplcr 
14, 15 and 16 (hereafter the cluster abbreviations Bo, Co 
and Tr will be used). These clusters contain a combined 
total of 64 O-type stars including one o f the most mas - 
sive and spectacular stars known, r\ Car l | FeinsteirJll995l) . 
Located at a distance of 2.2 kpc l|Txivmass!anlIl995|) . the 
Carina Nebula is an excellent region in which to study the 
effect massive stars have on their natal GMC. 
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The giant pillars are located in the relatively unstud- 
ied southern region of this nebula, at a greater distance 
from the most influential clusters. The stellar winds and 
radiation fields may be less destructive here, making the 
pillars prospective sites for ongoing star formation. It is 
the aim of this paper to investigate the nature of the in- 
terstellar medium within, and surrounding these pillars, 
and in particular to determine if there exists any evidence 
for ongoing star formation within them. 



2. The data 

2.1. MSX on-line database 

Images between 8-21 /im were obtained from the MSX 
database 1 , and cover four discrete bands; Band A (6.8- 
10.8 fmx), Band C (11.1-13.2 jum), Band D (13.5-15.9 jum) 
and Band E (18.2-25.1 /im). These bands will be hereafter 
referred to as the 8, 12, 14, and 21 /im bands respectively. 
The angular resolution of the images (FWHM) is 20" with 
an astrometric accuracy of ~ 2". 

Over the four MSX bands, the emission arises from 
either fluorescently excited polycyclic aromatic hydrocar- 
bon (PAH) molecules or black-body emission from heated 
dust grains (150-400K). Emission from the family of PAH 
features (6.2, 7.7, 8.7, 11.3, 12.7 and 16.4 /mi), as well as 
several emission plateaus caused by small PAH molecules 
(6-9 /im, 11-13 /im and 15-20 urn) are all contained within 
the 8, 12 and 14 /im bands IjVerstraete et al.ll200l| . The 
emission seen in the 21 /im band generally arises from 
heated dust, and is useful in revealing the locations of 
heavily obscured young stars. 

The MSX Point Source Catalogue (PSCvl.O) provided 
position and flux measurements for point sources within 
the Carina Nebula (over the extent shown in Fig. ^) . A 
total of 526 point sources were found, however, the major- 
ity are only detected in the higher sensitivity 8 /im band. 
For the discussion presented in Sect. 14.21 a subset of this 
larger catalogue was used. Sources were extracted from 
the PSC if they had reliable fluxes in each of the 8-, 14- 
and 21 /im bands (corresponding to quality Hags 3 or 4). 
This reduced the list to a total of 33 sources. 

2.2. MOST on-line database 

The radio continuum data presented here were obtained 
as part of the Molonglo Galactic Plane Survey (MGPS) 
made with the Molonglo Observatory Synthesis Telescope 
(MOST) 2 . The centre frequency of the MOST is 843 MHz, 
with a bandwidth of 3 MHz. The synthesized beam is 
43"x43"cosec|i $|. For more deta ils of the specifica tions of 
the MOST see lMillsl lll98lf) and lRobertsonl (Il99lh . 



1 For a full description of the MSX satellite, the astrophys- 
ical ex p eriments an d the o b serving tech ni ques see |MilLet_al 
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2.3. SEST observations 

2.3.1. Spectral line observations 

Spectral line observations were obtained toward four 
mid-IR condensations (G287 73-0.92, G287.8 4-0.82, 
G287.93-0.99 and G288.07-0.80: ISmitn et a,l1l2000h using 
the 15-m Swedish-ESO Submillimetre Telescope (SEST) 3 
at La Silla Observatory, Chile. The condensations were 
observed in six molecular lines: 12 CO(1-0), 12 CO(2-l), 
13 CO(2-l), C 18 O(l-0), CS(3-2) and CS(2-1). Details of 
the observations are given in Tabled 

Both the IRAM-built 230/115 GHz and the SESIS 
150/100 GHz receivers were used. Transitions were ob- 
served simultaneously where possible in single side band 
mode (SSB), with the receivers connected to a narrow- 
band (43 MHz) Acousto-Optical Spectrometer (AOS). 

Position-switching mode was used for all observations, 
with rj Car used as the signal-free reference position. An 
off-source position was taken after every second on-source 
integration. Chopper wheel calibration was performed ev- 
ery lOmins, to obtain atmosphere-corrected antenna tem- 
peratures. The telescope pointing and sub-reflector focus- 
ing were checked regularly using suitably bright nearby 
SiO masers. We estimate the pointing accuracy to be bet- 
ter than 10" and adopt the standard SEST value of 10 
per cent uncertainty in the temperature scale. 

Observations toward the condensations were obtained 
using a 20" pointing grid and intended to cover the peak 
of the emission. However, in several cases data were only 
obtained over the tip of much larger structures. 

All spectra were processed using GILDAS software 
(Bui sson et alJ Il999l) . Baselines were initially removed 
from each spectrum and the temperature scale converted 
to main-beam brightness temperature using the values for 
the main beam efficiencies (?/ m b) given in Tabled For com- 
parison between data sets with differing beam sizes, the 
data were smoothed to the larger beam-size on a similarly 
sized grid. The beam size of the SEST at each frequency, 
the average system temperatures during the observations, 
and the average rms noise per spectral channel are also 
given in Tabled 



2.3.2. Continuum observations 

Observations of 1.2-mm continuum emission were also 
obtained with the SEST using the 37 channel Imaging 
Bolometer Array (SIMBA). The FWHM of each element 
is 23" within the array, with a separation between each of 
44". The observations were conducted in July 2002 in the 
fast-scanning mode (80 arcsec s _1 ). Data for a single point 
source in the Carina Nebula (G287.87-l.36) were obtained 
using a map size of 480" x 240" (az x elev). To complete 
this map, 31 sub-scans were required with each sub-scan 
separated by 8". Each map took ~ 4 mins with the map re- 
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3 The SEST is operated jointly by ESO and the Swedish 
National Facility for Radio Astronomy, Onsala Space 
Observatory, Chalmers University of Technology. 
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Table 1. Observing parameters for the molecular line data obtained with the SEST. t- mt refers to the integration 
time, HPBW is the half-power beam-width, r] m b is the main beam efficiency, T S j, s is the average system temperature, 
and a is the average rms noise per spectral channel for the observations. 



Transition 


Frequency 


Date 


tint 


HPBW 


»?m6 


Tsys 


a 




GHz 




S 


arcsec 




K 


K 


1L! CO(1-0) 


115.217 


Oct 2000 


30 


45 


0.70 


529 


0.7 


12 CO(2-l) 


230.537 


Oct 2000 


30 


23 


0.50 


625 


0.7 


C 18 O(l-0) 


109.782 


Mar 2002 a 


120 


49 


0.72 


321 


0.4 


13 CO(2-l) 


220.398 


Mar 2002° 


120 


24 


0.52 


488 


0.5 


CS(2-1) 


97.981 


Jan 2002 


60 


52 


0.73 


337 


0.4 


CS(3-2) 


146.969 


Jan 2002 


60 


34 


0.66 


433 


0.5 



a The observations at these frequencies for source G287.93-0.99 were obtained in Oct 2000. 



peated 6 times to improve the signal to noise. The pointing 
and sub-reflector focusing were checked prior to the obser- 
vations with sky-dip calibrations performed both before 
and after the mapping. Observations of Uranus were also 
obtained for flux calibration. The average opacity during 
the observations was ~ 0.14. 

All maps were reduced by applying the opacity correc- 
tions, fitting and subtracting a baseline, and removing the 
correlated sky noise. Maps were flux calibrated using the 
conversion factor obtained from the Uranus observations. 
The rms noise level in the final map was 9 mJy. All data 
reduction was achieved within the MOPSI package 4 . 

2.4. 2MASS on-line database 

Near-IR images were obtained from the Two Micron All 
Sky Survey (2MASS) 5 Second Incremental Data Release. 
The images cover the three near-IR bands: J (1.24 /im), 
H (1.66 /im), and K s (2.16 /im) down to limiting magni- 
tudes of 15.8, 15.1 and 14.3 respectively. The resolution of 
the final images is 2". 

3. The nature of the pillars 

This section aims to identify, and determine the nature 
of, the pillars within the southern Carina Nebula. We will 
attempt to trace the extent of the gas and dust within 
the region, determine how the radiation from the nearby 
clusters is penetrating through the GMC, locate ioniza- 
tion fronts, and attempt to characterise the properties of 
regions where star formation may be occurring. 

3.1. Photodissociation Regions 

Fig. n shows the Carina Nebula in 8 /jm emission. This 
image reveals many regions of extended 8 /im emission in 
addition to several bright isolated features. This band is 
dominated by PAH emission and is tracing photodissoci- 
ation regions (PDRs). These occur within a visual extinc- 
tion (A v ) ~1 mag on the surface of molecular material, 

4 MOPSI was developed by R. Zylka (Astrophysikalisches 
Institut, Universitat zu Koln, Germany). 

5 see http://www.ipac.caltech.edu/2mass/index.html 



where th e UV radiation has penetrate d and is heating the 
material l|Hollenbach fc Tielenslll999j) . 

Evident to the northwest of this nebula is a prominent 
8 /im feature which corresponds to a component of the 
GMC known as the northern molecular cloud. The bright 
emission south of r\ Car traces the edge of the southern 
molecular cloud. Locat ed further to the s outh is the gi- 
ant pillar identified bv lSmith et al These clouds 
are al l identified in the la rge scale CO(4-3) and [CI] sur- 
vey of lZhang et alJ (|200lf> . These transitions are tracers of 
heated gas (~50K), dense gas (n ~10 5 cm -3 ), and PDRs, 
respectively, and confirm that the 8 /jm emission is tracing 
the outer PDR layers of the molecular clouds. 

Fig. 1H also illus trates the mid-IR sources identi- 
fied by ISmith et al] (2000; G287.73-0.92, G287.84-0.82, 
G287.93-0.99 and G288.07-0.80). Many other point 
sources are also found across the complex, and are located 
at the tips of more extended 8 /im emission (e.g. G287.57- 
0.91 and G287.88-0.93). Their morphology is consistent 
with a dense core being irradiated by Tr 14 and Tr 16. 

3.2. Visual extinction 

Extinction maps were constructed for the southern Carina 
Nebula using stars within the 2MASS PSC. Estimates of 
the visual extinction were obtained for each star detected 
in the JHK s -bands, by dereddening each to an average 
main sequence colour in the (J-H) vs. (H-K s ) plane. The 
visual extinctions obtained were then used to construct a 
global extinction map for the region (contours on Fig. ^ 
the extent of the map is marked with the dotted line) . 

Clearly seen as a region of high visual extinction is 
the giant pillar, where the morphology of the extinction 
matches the 8 /im emission extremely well. In addition to 
this giant pillar, other regions of high visual extinction are 
also found to correspond to 8 /im features. 

The visual extinction across the giant pillar ranges 
from 7-15 mag and can be used to estimate the 
H2 column density within it. Applying the rela tion, 
N( 13 CO)~ 2.7xl0 15 A v fcm" 2 : iFrerking et al.ll98l and 
convertin g to N(H?) using t he constant of proportionality 
of 4xl0 5 IjLada et al.ll"l994h . the H 2 column density was 
found to range between 9-19xl0 21 cm~ 2 . 
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Fig. 1. The Carina Nebula in emission (logarithmic grey-scale) with contours of visual extinction (levels are 

7.5, 9.5, 11.5, 13.5 and 15.5 mag). Marked on the figure are the locations of r] Car (+), Tr 14 (x), the southern and 
northern molecular clouds, and the additional giant pillar and molecular cloud identified by Zhang et al. (2001). The 
dotted region outlines the extent of the extinction map. The four condensations labelled (in white) G287.73-0.92, 
G287.84-0.82, G287.93-0.99 and G288.07-0.80 are discussed in Sect.ESl The point sources labelled in black mark the 
sources discussed in Sect. 0] 



3.3. Heated dust 

Fig. shows the 21 /xm and 843 MHz emission across the 
Carina Nebula. Of note, is the striking morphological sim- 
ilarities of these emission features across the complex. 

Strong 21 fim emission surrounds rj Car and Tr 14 
with diffuse emission found in the south. Several 
bright point sources are also evident, in particular the 
bright source associated with the 8 fim condensation 
G287.84-0.82. As mentioned previously, extended emis- 
sion seen within the 21 /im band generally arises from 
heated dust, while 21 /im point sources identify deeply 



embedded stars. By matching the expected surface bright- 
ness of an optically thick black-body with the observed 
flux from regions showing strong 21 fim emission, we find 
a dust temperature of 42-44 K. This temperature is in 
good agreement with other esti mates (40K) obtained from 
IRAS data in the Tr 14 region jHarvev et alJll979|) . 

3.4. Ionization fronts 

Radio continuum emission is widespread throughout the 
Carina Nebula, with the brightest concentrations known 
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Fig. 2. The Carina Nebula in 21 fim emission (logarithmic grey-scale) with contours of 843 MHz radio continuum 
emission (levels are 0.04, 0.14, 0.24, 0.34, 0.5, 0.7, 0.9, 1.1, 1.3, 1.5 Jy/beam). Marked on this figure is the location 
of r\ Car as well as several ionization fronts (including Car I and Car II), and the four condensations; G287.73-0.92, 
G287.84-0.82, G287.93-0.99 and G288.07-0.80. 



as Car I and Car TT <TBr"ooks et a.lJl20mt IWhiteoaklll 994t 
lRetallacklll983lGardner & Morimotolll968|) . 

Extending from Car I is a long filament of weaker emis- 
sion. Similar features are also seen in the southern part of 
the nebula, and are labelled on Fig. [5] All of these features 
are coincident with 21 /an emission, and in the majority of 
cases, are found adjacent to regions of bright 8 fim emis- 
sion and peaks in the visual extinction map. This mor- 
phology is consistent with the ionization fronts heating 
the dust and creating the widespread PDRs. 

Compact sources are also evident across the com- 
plex. The majority of these have no mid-IR counter- 
parts (marked as 'point sources' on Fig. |2J), suggest- 



ing they correspond t o unresolved background galaxies 
l|Cohen k Greenl l200l[) . The exception is G287.84-0.82, 
which has an extremely strong mid-IR counterpart. This 
object is likely to be a compact Hn region. 

3.5. Properties of the four bright mid-IR condensations 

Molecular material was detected toward each of the 
four mid-IR condensations, G287.73-0.92, G287.84-0.82, 
G287.93-0.99 and G288.07-0.80. Their integrated 
12 CO(2-l) intensity maps, show the morphology of the 
molecular structures match well with the corresponding 
8 /im emission (Fig. EJ), further confirming the presence 
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Fig. 3. 8 /im emission in grey-scale with contours of integrated 12 CO(2~l) emission for the condensations, (a) 
G287.73-0.92: Integrated from -25 to -32 kms" 1 , levels are 20, 35, 50, 65, 80, gSKkms" 1 ; (b) G287.84-0.82: Integrated 
from -12 to -16 kms" 1 , levels are 20, 40, 60, 80, 100, 140, lSOKkms" 1 ; (c) G287.93-0.99: Integrated from -18 to - 
22 kms -1 , levels are 30, 50, 70, 90, 110, lSOKkms" 1 ; (d) G288.07-0.80: Integrated from -19 to -24 kms" 1 , levels are 
15, 30, 45, 60, 75, 90, lOSKkms" 1 . 



Table 2. Properties of the spatially integrated profiles seen toward the condensations. These were obtained from 
Gaussian fits to the profiles shown in Fig.0] The parameters include; the central velocity (V, km s -1 ), peak temperature 
(T, K) and line-width (AV, kms" 1 ). 
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Table 3. Physical parameters for the condensations derived from the properties given in Table |2 Listed are: the 
excitation temperature ( 12 T 0X ) measured from the 12 CO(2-l) peak temperature; the ratio of 12 CO(2-l) and (1-0) 
integrated intensities ( 12 R); the ratio of 13 CO(2-l) and 12 CO(2-l) integrated intensities ( 12 ' 13 R); the ratio of CS(3-2) 
and CS(2-1) integrated intensities ( CS R); 13 CO(2-l) optical depth ( 13 r); 13 CO(2-l) column density (N( 13 CO)); H 2 
column density (N(H2)); mass estimated under LTE conditions (Mltb); Virial mass estimate (M V i r i a i) and the average 
H 2 density (n(H 2 )). 
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of PDRs on their surfaces. For G288.07-0.80 however, the 
geometry of the molecular material differs considerably 
from the 8/im emission. While the morphology of the 
molecular gas is consistent with irradiation from the 
direction of the clusters Tr 14 and Tr 16, the PDR 
appears to point in a different direction. Interestingly, 
this corresponds to the direction of the Bo 11 cluster. 

Fig. 0] shows the line profiles for each condensation, 
in each of the molecular transitions observed. The peak 
temperature (T) , central velocity (V) and line- width ( AV) 



for these profiles were determined from Gaussian fits, the 
results of which are given in Table [21 

The molecular condensations show a range in veloc- 
ities of ~15kms _1 . This is consistent with the velocity 
dispersion of molecular clumps, and motion of the ionized 
gas, observed within t he Keyhole Nebula (~20kms , 
ICox k Bronfman! Il995|) . This suggests there is also a 
strong interaction between the winds and radiation from 
the stellar clusters and these molecular condensations. 
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Fig. 4. Line profiles spatially integrated over each con- 
densation. Left to right: G287.73-0.92, G287.84-0.82, 
G287.93-0.99 and G288.07-0.80. Top to bottom: 
12 CO(2-l) and 12 CO(1-0); 13 CO(2-l) and C 18 0(1- 
0); CS(3-2) and CS(2-1) emission. The C 18 O(l-0) 
profiles are multiplied by 5 to make them visible. Note 
the different temperature scales and velocity ranges. 



Estimates of the physical properties of the conden- 
sations were determined from the profiles under the as- 
sumption of local thermal equilibrium (LTE). Values for 
the excitation temperatures, line ratios, opacities, masses, 
and average hydrogen densities for each condensations are 
listed in Table |31 The ratio of the integrated intensity of 
the 12 CO(2-l) and (1-0) emission ranged between 1-2, 
and is consistent with optically thick 12 CO emission. In 
addition, excitation temperatures were found to range be- 
tween 20-40 K. 

Ratios of the 12 CO(2-l) and 13 CO(2-l) emission 
ranged between 2-5, suggesting that the 13 CO(2-l) tran- 
sition is optically thin. To estimate the opacity and column 
density of this material, a constant excitation temperature 
of 35 K was used. Opacities ranged from 0.16-0.4, confirm- 
ing the emission is optically thin. H2 column densities 
were determined to be 7-22 xlO 21 cm -2 , and are consis- 
tent with values determined when using the estimated A v 
derived from the star colours (9-19xl0 21 cm -2 ). 

Estimates of the masses were calculated using the col- 
umn densities determined from the 13 CO emission (LTE 
mass). This calculation assumes the integrated intensity of 
the emission is proportional to the H2 c olumn density wit h 
a constant of proportionality of 4x 10 5 llLada et alJIllffli . 
Masses ranged from 10-45 M Q for the condensations, how- 
ever, these may be under-estimated by as much as a factor 
of 3. The Virial mass was estimated using the line- width 
of the 13 CO line, which is assumed to be a measure of 
the overall motion of the gas and hence the mass. Values 
of the Virial mass for the condensations ranged between 
80-200 Mq. These are significantly greater than the LTE 
masses (even considering the latter may be underestimat- 
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Fig. 5. K s -band image for the southern Carina Nebula 
with contours representing the stellar density distribution 
for stars with K corr < 12 mag (thick black: levels are 2.6, 
3, 3.4, 3.8 stars/arcmin 2 ) and the visual extinction map 
(thin black: levels are 7.5, 9.5, 11.5, 13.5 and 15.5 mag). 
The cluster Bo 11 is marked, along with the condensations 
G287.84-0.82 and G287.93-0.99. 



ing the true masses), implying that the condensations are 
probably not gravitationally bound but instead confined 
by external pressure. 

The average H 2 density was determined assuming the 
size along the line-of-sight was the average of the extent in 
right ascension and declination. Estimates ranged between 
4-17 x 10 3 cm~ 3 and represent lower limits to the ac- 
tual values. In addition, the presence of CS emission sug- 
gests the material has gas densities in excess of 10 4 cm -3 . 
The properties described here are consistent with molecu- 
lar cores found in other regions of massive star formation 
(Zinchenko et al 1998). 

4. Star formation within the pillars 

This section investigates the star formation activity across 
the Carina Nebula, by indentifying candidate embedded 
clusters and young stellar objects (YSOs). 



4.1. Young clusters 

All stars within the 2MASS PSC were dereddened. As pre- 
viously discussed, those detected in the JH and K s -bands 



8 



Rathborne et al.: The giant pillars of the Carina Nebula 



2 



0.5 



_ 1 1 ' 


,u 


- 

MYSOs? 












■ _ 

4 A 

■ 

j- 200K 


A 




1 


A 




* J * 




-^A 






AA 






A 


2 1 A 


' 5000K A 


A 






▲ 

A Compact HI 
regions? 



^21 

Fig. 6. Colour-colour plot for sources in the Carina 
Nebula. Included on this are the location of a black- 
body (solid line) and the limits derived by Lumsden et 
al. (2002) for MYSOs and compact Hn regions (dashed 
lines). Plotted here are: point sources with MSX data 
(filled triangles), point sources with MSX and IRAS de- 
tections (open triangles), MSX sources associated with 
diffuse emission (filled squares) and the source identified 
by Megeath et al. (1996) (open asterisk) and Rathborne et 
al. (2002) (filled asterisk). The numbers mark the sources 
discussed in Sect. 14.31 



were dereddened using an average main sequence colour. 
For those stars only detected in the H- and/or K s -bands, 
the extinction was estimated using nearby stars in the 
global extinction map. From these dereddened K s -band 
magnitudes (K corr ) a stellar number density map was cre- 
ated for stars with K corr >12 mag. Choosing this cut-off 
corresponds to the completeness limit from the 2MASS 
photometry for an A v of 25 mag. 

With the foreground extinction removed, the brightest 
sources will correspond to the most massive and youngest 
objects. Young sources display a near-IR excess, making 
them brighter in the K s -band than main-sequence stars 
whose emission is purely from a stellar photosphere. As 
stars tend to form in clusters, any over-densities within the 
stellar density distribution will potentially reveal the lo- 
cation of clusters. Fig. displays, as contours, the derived 
stellar density, and visual extinction maps for the south- 
ern Carina Nebula. Many over-densities are seen here and 
are candidate young clusters. 

Coincident with a peak in the stellar density distribu- 
tion is the known cluster Bo 11 and interestingly, two of 
the molecular condensations discussed in Sect. 13.51 The 
source G287.84-0.82 reveals, in 2MASS images, extended 
IR emission and a tightly packed group of sources, many 
of which are only visible in the K s -band image (this source 
will be discussed further in later sections). 

Many other candidate clusters are identified from 
this analysis, but display no obvious clusters within the 



2MASS images. Many are associated with peaks in the vi- 
sual extinction maps, therefore making the identification 
of a cluster from the images difficult. While clusters are 
located across the region, it is those associated with the 
visual extinction peaks that are potentially the most inter- 
esting. This correspondence suggests the clusters may be 
embedded within the pillars. Additional data is required, 
however, to determine the exact nature of many of these 
candidate clusters. 

4.2. Identification of candidate YSOs 

Recently, iLumsden et alJ (j2002^l derived criteria for iden- 
tifying various galactic plane sources from fluxes obtained 
within the MSX bands. The main focus of their study was 
to separate massive young stellar objects (MYSOs) from 
sources with similar IR colours, in particular, compact H n 
regions, young compact planetary nebulae, and very dusty 
evolved stars. From mid-IR colour-colour diagrams, they 
find that young objects satisfy F2i/Fg > 2, with MYSOs 
having F14/F12 > 1, while compact Hn regions have 
F14/F12 < 1. These limits were determined from the loca- 
tion of a known sample of these objects. The sample of H 11 
regions originate d from sources identified within radio con- 
tinuum surveys llWood fc Churchwelll Il989t iKurtz et alJ 
119941 IWalsh et a,l.lll998h. For the MYSOs, the list was 
compiled by Lumsden et a ll2002l) . and consisted of ob- 
jects with strong emission in the longer MSX bands, and 
with featureless mid-IR spectra, with the possible excep- 
tion of silicate absorption. In this section, we use these cri- 
teria to investigate the nature of the point sources across 
the Carina Nebula in order to determine if, and where, 
any MYSOs and compact Hn regions may exist. 

Fig. shows the colour-colour diagram for all sources 
with reliable detections in the MSX and/or IRAS PSCs. 
Also included on this plot is the position o f a bl ack- 
body and the limits derived bv lLumsden et aD |2002) for 
MYSOs and compact H 11 regions. 

In this colour-colour plane, sources with F14/F12 < 1 
and F21/F8 < 2 correspond to evolved stars. The remain- 
ing sources that satisfy F2i/Fg > 2, are candidate young 
objects, 13 of which fall within the compact H 11 region 
while 9 are potentially MYSOs. Visual inspection of the 
images reveal the sources displaying an F2i/Fg > 10 are 
associated with emission and are not point sources (these 
are marked as filled squares in Fig. [SJ). Thus, their identi- 
fication as MYSO objects is dubious. Two MYSO candi- 
dates remain and are located in the south of the nebula 
(labelled on Fig. □ as G287.87-l.36 and G287.73-l.01). 
These appear isolated from any extended mid-IR, radio 
continuum emission, or the molecular condensations. A 
candidate compact H 11 region is, however, associated with 
one of the molecular condensations (G287.84-0.82). 

Two sources identified within this analysis have pre- 
viously been identified as YSOs and include a bright- 
rimmed globule at the edge of the southern molecu- 
lar cloud (marked with an open asterisk; G287.63-0.72; 
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Fig. 7. Spectral energy distributions for the candidate YSOs; G287.52-0.41 (source 1), G287.84-0.82 (sources 2), 
G287.88-0.93 (source 3) and G287.87-l.36 (source 4). Open squares represent 2MASS fluxes, stars represent MSX 
fluxes, while circles mark the IRAS fluxes. The solid line in each plot is a fit to these points using a modified two 
component black-body function (see text for details). The 8^m flux is excluded from the fit for G287.84-0.82 as it is 
known to be dominated by emission from PAH molecules. The near-IR points were not included in the fits as these 
likely represent emission from an extincted hot component rather than dust emission. 



Table 4. The derived parameters of the two-component grey-body fits to the SEDs for G287.52-0.41, G287.84-0.82 
G287.88-0.93 and G287.87-l.36 as shown in Fig. The parameters correspond to the temperatures of the outer 
(T ou ter) and inner (Ti nner ) dust components, the inner (Rinner) and outer (Router) radii of the dust shell and estimates 
of the gas mass (M gas ), Luminosity derived from these and the corresponding spectral types of the exciting star. The 
numbers in brackets next to the source names refer to the source label from Fig. 



Source 


Touter 


Tinner 


Runner 


-Router 


M gas a 


Luminosity 


Spectral 




K 


K 


AU 


AU 


M Q 


10 4 L 


Type 


G287.52-0.41 (1) 


80 


320 


10 


850 


16 


5.1 


08.5 


G287.84-0.82 (2) 


60 


300 


57 


7200 


1200 


370 


cluster? 


G287.88-0.93 (3) 


90 


330 


14 


570 


7.3 


3.4 


B0 


G287.87-l.36 (4) 


80 


360 


14 


960 


12 


3.3 


B0 



a The mass of G287.87-l.36 was determined from the observed flux at 1.2-mm. Estimates of the 
1.2-mm flux for the other sources were obtained from the grey-body fit. 



iMeeeath et al.ll99fj) and an embedded 09- star (filled as- 
terisk; G287.47-0.54; lRathborne et al.l2002f) . 

4.3. Continuum SEDs 

To investigate the nature of the YSO candidates, spec- 
tral energy distributions (SEDs) were constructed. Of the 
14 candidates, only 4 of these had reliable and coincident 
detections within both the MSX and IRAS PSCs. The 
first of these sources is located in the northern part of the 
nebula (G287.52-0.41), the second is the source associ- 
ated with the molecular condensation (G287.84-0.82), the 
third is located at the tip of a giant pillar (G287.93-0.99), 
while the remaining source is located in the south, and is 
a MYSO candidate (G287.87-l.36). These are discussed in 
more detail in Sect. 14.51 

Fig. shows the continuum SEDs for these sources. 
Included on these plots are the fluxes obtained from the 
IRAS PSC, as well as those estimated from the MSX 
and 2MASS images. An additional 1.2 mm flux is in- 



cluded for the source G287.87-l.36. The solid line in each 
plot is a fit to the data using a grey-body of the form 
B J/ (T)(l-exp(-r )y ))n, where B y (T) is the Planck function 
at a temperature T, t v is the dust optical depth, and £1 
is the solid angle subtended by the dust shell. The opac- 
ity was assumed to have the form r = (yjv^ where v 
is the frequency at which the emission becomes optically 
thic k, with set to 2 (determ ined from tabulated values 
of lOssenkonf fc Hennindll994 - A two-component model 
was fit to the data, because a single temperature com- 
ponent was insufficient to describe the observed emission. 
Parameters derived from these fits, including the temper- 
atures and radii of the two dust components, mass, lumi- 
nosity, and spectral type of the central star, are listed in 
Table H 

Observations at a frequency in which the emission is 
optically thin allows a determination of the mass of the 
material. For the source G287.87-l.36 additional data at 
1.2-mm were obtained which, when included in fit, re- 
vealed a vq of 4.0xlO n Hz (0.75mm). Thus, the emission 
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at 1.2 mm is optically thin and was therefore used to es- 
timate the mass. The expression M gas = fld f"g ( Td ) was 
used, where M gas is the total mass of an isothermal dust 
source, Sj, is the observed flux density at an optically thin 
frequency v, n v is the mass absorption coefficient, Kd g is 
the dust-to-gas mass ratio and B,ATd) is the Plan ck func- 
tion at the dust temperature (jChini et al.l Il987^ . Using 
values of R d„ ~ 0.01, a dust opacity of 0.1 m 2 /kg (deter- 
mined from lOssenkopf fc Henningl ll994). and assuming a 
distance of 2.2 kpc, the resulting mass was found to be 
12 Mq. The mass estimates represent a lower limit, as 
the dust opacity, and temperatures are uncertain, and are 
most likely lower than quoted here. 

Luminosity estimates suggest the exciting stars for 
these sources are OB-stars. The high mass and luminosity 
derived for G287.84-0.82 suggest this object may harbour 
a cluster of stars. 

4.4. Near-IR three-colour images 

Fig.[H]shows near-IR three-colour images constructed from 
the 2MASS data, for the 14 YSO candidates plus one 
source of interest (G287.86-0.82). This source was selected 
because of its extremely red near-IR colour, similar to the 
MYSO candidate G287.87-l.36. These images reveal any 
near-IR counterparts to the mid-IR sources, in particular, 
note the prominant cluster of reddened stars associated 
with G287.84-0.82. 

Near-IR images of the sources G287.73-l.01 and 
G287.80-0.56, reveal bright stellar objects, the former be- 
ing a ssociated with a known 07-star (HD 93222; W albornl 
Il972f) . Although these objects lie in the regime for candi- 
date MYSOs and compact Hu regions respectively, they 
most likely represent older stars, as they are clearly seen 
in the near-IR and, for the object G287.73-l.01, also at 
visible wavelengths. As a result, these are not considered 
to be candidate YSOs. 

Several sources identified as compact Hu regions, are 
associated with diffuse emission in the near-IR images. 
These objects (G287.51-0.49, G287.22-0.53, G287.47-0.54, 
G287.63-0.72, G287.75-0.66, G287.76-0.87 and G287.68- 
1.09) all show extended K s -band emission, with several 
reddened sources located at the peak of the 8 /jm emis- 
sion. In addition, all seven of these objects have extended 
843 MHz radio continuum emission associated with them. 

Four of these objects (G287.51-0.49, G287.22-0.53, 
G287.47-0.54, and G287.63-0.72) contain multiple sources 
and are possibly embedded clusters. These objects were 
all located external to the stellar density map of 
Sect.O The sources G287.4 7-0.54 and G287. 63-0. 72 have 
been discussed previo usly in lRathborne et all l|2002(l and 
iMeeeath et"aT1 (|l996^l . 

Near-IR images of the remaining three sources, 
G287.75-0.66, G287.76-0.87 and G287.68-l.09, show no 
obvious embedded stars, with the stellar colours most 
likely a result of interstellar extinction. 
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Fig. 9. The IRAS LRS spectrum of G287.84-0.82 (open 
squares) with error bars, and comparisons with photome- 
try from the IRAS PSC (filled triangles) and MSX PSC 
(filled squares) and from the MSX images (open circles). 

4.5. Synopses of the most interesting YSO candidates 

4.5.1. G287.52-0.41 

This source is located in the north of the Carina Nebula 
and has mid-IR colours indicative of a compact H n region 
(source 1 in Fig. |5J). Consistent with this, is the detection 
of 843 MHz radio continuum emission toward this object. 
Mass and luminosity estimates suggests the presence of an 
embedded 08.5-star. Furthermore, in the near-IR, a single 
reddened object is seen and is coincident with the MSX 
and IRAS point sources. 

4.5.2. The cluster within G287.84-0.82 

This is the only case within the current dataset, for which 
a clearly visible cluster of young stars has been found em- 
bedded within a molecular condensatio n (this cluster was 
first identified bv iDutra fc Bicall200lj) . It is the bright- 
est point source in the mid-IR and 843 MHz radio contin- 
uum maps within the southern Carina Nebula. The near- 
IR three colour image for this cluster, reveals many em- 
bedded stars, with several detected only at the K s -band. 
These sources represent the most interesting and youngest 
within the cluster. 

In the mid-IR colour-colour plane, this source (num- 
ber 2), lies within the region for compact Hu regions. As 
expected, it shows strong, compact radio continuum emis- 
sion. The mass and luminosity derived for this source, sug- 
gests it contains at least one massive star within a cluster 
of lower mass stars. This is consistent with previous stud- 
ies which su ggest a Bl star and reflection nebula exist here 
(ALS 1883: iHerbstl Il975h . However, 6.7 GHz methanol 
maser emis sion was not detected toward this object in 
a survey bv lMacLeod et alJ l|1998|) . 

An averaged IRAS Low Resolution Spectrometer 
(LRS) spectra was also obtained for this object (Fig. EJ). 
Within the sizeable LRS aperture (5' at 8-13 fim and 



Rathborne et al.: The giant pillars of the Carina Nebula 



11 




Fig. 8. A collection of near-IR three-colour images for candidate YSOs within the Carina Nebula (J-band in blue, 
H-band in green and K s -band in red). The images are presented in the order in which they are discussed in Sect. 14.51 
Included also are yellow ellipses marking the location of IRAS point sources and contours of 8 fim emission. Levels 
are, G287.52-0.41: 0.5, 1; G287.84-0.82: 1, 2, 4, 8, 10; G287.88-0.93: 0.5, 1; G287.57-0.91: 0.5, 1, 1.5; G287.87-l.36: 0.5, 
1, 2; G287.86-0.82: 0.5, 0.7; G287.73-l.01: 0.5, 1, 2; G287.80-0.56: 1, 2; G287.51-0.49: 1, 2, 3; G287.22-0.53: 0.8, 1.3; 
G287.47-0.54: 2, 4, 6; G287.63-0.72: 1, 2, 3; G287.75-0.66: 2, 3, 4; G287.76-0.87: 1, 2; G287.68-l.09: 1, 1.5. All are in 
units of 10~ 5 W m~ 2 sr _1 . 
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Fig. 10. The IRAS LRS spectrum of G287.87-l.36 (open 
squares) with error bars, and comparisons with photom- 
etry from the IRAS (filled triangles) and MSX (filled 
squares) PSCs and from the MSX images (open circles). 

7.5' at 11-23 fim) strong PAH emission is seen at 7.7 and 
11.3 /im, confirming our previous identification of PAHs 
on the basis of the MSX and 12 CO morphology. Flux 
densities estimated from the MSX images and PSC, and 
IRAS PSC, are also included on this plot. The obvious dis- 
crepancies between flux densities derived from the LRS, 
the MSX image, and those from the MSX PSC, reflect 
the complexity of the region, with its widespread mid-IR 
emission, and the difficulties in choosing an appropriate 
background level to subtract from the dominant mid-IR 
peak. 

4.5.3. G287.88-0.93 and G287.57-0.91 

The mid-IR morphology of these two sources show fila- 
mentary features with point sources at their tips. They 
were included within the colour-colour plot of Fig. and 
fall within the regime for compact Hn regions. However, 
no strong 843 MHz radio continuum emission is associated 
with them. The source G287.88-0.93 does, however, con- 
tain an IRAS point source, with its SED corresponding to 
an embedded BO star. A star of this spectral type would 
produce a compact H n region, albeit weak, and possibly 
below the detection limits of the MOST considering the 
distance to the nebula. 

The near-IR image for G287.88-0.93 shows a group of 
three stars falling within the peak of the 8 /im emission. 
Near-IR images of G287.57-0.91, reveal a bright K s -band 
source located at the tip of the filamentary 8 /im feature. 

4.5.4. G287.87-l.36 

This source is the most likely to correspond to a bona-fide 
MYSO, as it lies within this regime in the mid-IR colour- 
colour diagram, and has a corresponding IRAS point 
source. The luminosity derived from its SED suggests an 
embedded BO-star. An extremely red object exists in the 



near-IR three colour image, and is coincident with an un- 
resolved mid-IR and 1.2-mm continuum source. The IRAS 
LRS for this source (Fig. llOjl . clearly shows a deep silicate 
absorption profile at 10 /im indicative of high extinction. 

In addition, iBronfman et all l)l996|) detected CS (2- 
1) line emission toward this object in their survey 
of IRAS sources having colours of compact Hn re- 
gions. It appears in their table 1 as IRAS 10423-6011, 
with V isr =13.7kms- 1 , T m6 =0.84K and AV=2.7kms- 1 . 
Because the central velocity is blue-shifted by more than 
15kms _1 with respect to the other molecular condensa- 
tions here, it is not clear if this object is actually associated 
with the nebula. 

4.5.5. G287.86-0.82 

Located to the east of the cluster G287.84-0.82 is the ob- 
ject G287.86-0.82. This source shows weak 8 /im emission, 
with a strong K s -band source. Located nearby, is a second 
object whose emission peaks in the H-band with extended 
K s -band emission (green in the near-IR three colour im- 
age). Such emission either arises from a reflection nebula 
or shocked gas. Detailed spectra are needed to determine 
the exact nature of this object. 

4.6. Ongoing star formation within the pillars 

The sources identified here significantly increase the num- 
ber of young objects found within the Carina Nebula. 
Previous studies located 6 embedded OB-stars in the 
northern molecular cloud complex associa ted with Tr 14 
l|Brooks et, al.ll200lURathhorne et all2002j >. while a single 
BO star is seen at the edge of the southern molecular cloud 
l|Megeath et al.lll996h . 

To investigate the age distribution of the star forma- 
tion within this GMC, the recent rate of massive star 
formation (ip) has been estimated. Studies of the stel- 
lar population within the clusters Tr 14 and Tr 16 find 
that the massive stars formed within the past 1 -3 Myr 
l|DeGioia-Eastwood et al1l200li iTapia et al.ll2003j) . These 
clusters have a population of ~ 33 O-type stars, which 
implies a ip of ~ 1-3 xlO -5 O-stars yr _1 . 

Mass and luminosity estimates exist for ten young ob- 
jects within the Carina Nebula (the seven previously iden- 
tified and the three additional sources identified here) . Age 
estimates were calculations from their Kelvin-Helmholtz 
time, assuming they consist of a single star. Ages were 
found to range between 1-5 xlO 5 yrs, which for these 
10 sources, suggests tp to be ~ 3xl0 -5 O-stars yr" 1 . 
This value is consistent with the rate determined from the 
known optical O-stars and is expected if the rate of star 
formation has been constant over the age of the nebula. 
While this is possible, these two epochs may instead repre- 
sent two isolated bursts of star formation. Distinguishing 
between these two mechanisms is difficult with the limited 
data presented here. 
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The rate of star formation for the Carina Nebula is 
comparable with estimates for two other well studied re- 
gions of massive star formation. For instance, NGC 3606 
contains ~ 50 O-t ype stars which have fo rmed over the 
last 2-3 Myrs fsee lEisenhauer et al .11 19981 and references 
therein). This suggests a ip of ~ 1-3 xl0~ 5 O-stars yr" 1 . 
While for 30 Dor, at a similar age, and containing ~ 100 
O-type stars IWalborn fc Blades! l|l997|) . the rate is esti- 
mated to be ~ 3-5 xl0~ 5 O-stars yr -1 . 



5. Conclusions 

We have undertaken a multi-wavelength study incorpo- 
rating data from 2MASS, MSX, IRAS, MOST, and the 
SEST, to investigate the nature of, and search for star 
formation within, the giant pillars of the Carina Nebula. 
Our main results and conclusions are summarised below. 

5.1. Molecular clouds, PDRs and ionization fronts 

Emission from the 8 /im MSX band outlines the known 
molecular clouds and giant pillars within the Carina 
Nebula, and pinpoints regions where the UV radiation is 
penetrating the molecular material and forming PDRs. 
Visual extinction maps match extremely well the 8 /im 
emission tracing the dense gas associated with the pillars. 

Interestingly, the 21 /im and 843 MHz radio continuum 
emission match extremely well and are located along the 
edges of the PDRs closest to the most influential clusters. 
Emission within these bands reveals heated ~40 K dust 
and ionization fronts. The geometry of the largest pillar is 
consistent with interactions from the nearby massive stars 
carving the molecular material around a dense core. 

Bright condensations located at the tips of the giant 
pillars, were found to be externally heated, with PDRs 
located along the edges in the direction of the massive 
clusters. The properties of the molecular material suggest 
they have the potential to be massive star forming cores. 

5.2. Evidence for recent star formation activity 

To search for evidence of star formation activity a stel- 
lar number density map was used. This was derived from 
the dereddened K s -band magnitudes for sources in the 
2MASS PSC. Many candidate young clusters were iden- 
tified within this map, several of which appear to be re- 
lated to, and possibly embedded within, the giant pillars. 
Interestingly, the brightest IR and radio continuum source 
in the region contained the only case of a clearly visible 
cluster within the 2MASS images. 

A candidate MYSO and several compact H n regions 
were identified across the nebula using mid-IR colour crite- 
ria. SEDs were constructed to study these objects further, 
and reveal the exciting stars in several cases correspond to 
OB-stars. TableQOlists a summary of all the candidates. 



5.3. Triggered star formation in the pillars? 

The results presented here clearly show that the large 
molecular clouds within the Carina Nebula are being 
strongly affected by the intense stellar winds and harsh 
radiation fields from the nearby massive star clusters. For 
instance, giant pillars are forming as a result of these in- 
teractions, as dense cores shield the surrounding material. 
What still remains unclear however, is the effect these in- 
teractions have on the ongoing star formation within this 
GMC. Young clusters and candidate MYSOs are found 
both within, and external to, the giant pillars. In addi- 
tion, while star formating sites are found scattered across 
the region, there are several locations where the intense 
interactions from the known clusters are potentially trig- 
gering the activity. It is also likely that the rate of star 
formation has been constant within the nebula since the 
birth of the massive star clusters. 
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Appendix 

Table [5] lists the coordinates and summarizes the proper- 
ties of the sources discussed within this work. 
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